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Abstract 

Based upon the formalism recently developed by one of us (MS), we analytically 
perform the post-Newtonian expansion of gravitational waves from a test particle in 
circular orbit of radius tq around a Schwarzschild black hole of mass M. We calculate 
gravitational wave forms and luminosity up to t> 8 order beyond Newtonian, where 
v = (M/ro) 1//2 . In particular, we give the exact analytical values of the coefficients 
of In v terms at v 6 and v 8 orders in the luminosity and confirm the numerical values 
obtained previously by the other of us (HT) and Nakamura. Our result is valid in the 
small mass limit of one body and gives an important guideline for the gravitational 
wave physics of coalescing compact binaries. 



1 Introduction 



Among the possible sources of gravitational waves, coalescing compact binaries are the 
most promising candidates which can be detected by the near-future laser interferometric 
gravitational wave detectors such as LIGO [||] and VIRGO 0. One reason is that we 
expect such events to occur 3/yr within 200Mpc |§. The other reason is that we expect 
enough amplitude of gravitational waves to be detected by LIGO and VIRGO if such 
events occur. 

Gravitational radiation from coalescing compact binaries contains rich information 
about physics of neutron stars, cosmological parameters, a test of general relativity and 
so on. Such information can be extracted out from the gravitational wave form by the 
matched filtering technique, that is, by cross-correlating the incoming noisy signal with 
theoretical templates. If the signal and the templates get out of phase with each other 
by one cycle as the waves sweep through the LIGO/ VIRGO band, their cross correlation 
will be significantly reduced. This means that it is important to construct theoretical 
templates which are accurate to better than one cycle during entire sweep through the 
LIGO/ VIRGO band 0. Thus, much effort has been recently made to construct accurate 
theoretical templates || 0, |, (|- 

To construct theoretical templates, the post-Newtonian approximations are usually 
employed to solve the Einstein equations. However, based on numerical calculations of the 
gravitational radiation from a particle in circular orbit around a non-rotating black hole, 
Cutler et al. [|TU[ showed that evaluation of the gravitational wave luminosity to a post- 



Newtonian order much higher than presently achieved level will be required to construct 
the templates. Then in order to find out the necessary post-Newtonian order, the same 



problem was investigated by Tagoshi and Nakamura ||TT| with much higher accuracy. 
They calculated the coefficients of the post-Newtonian expansion of the gravitational 
wave luminosity to (post) 4 -Newtonian order (i.e.,0(u 8 ) beyond Newtonian) and concluded 
that the accuracy to at least (post) 3 -Newtonian order is required for the construction of 
theoretical templates. In addition, they found logarithmic terms in the luminosity at 
(post) 3 and (post) 4 -Newtonian orders. 

It is then highly desirable to reproduce these results in a purely analytical way, that 
is, to derive the exact analytical expressions for the coefficients of the post-Newtonian 
expansion. Poisson |TSj first developed such a method and calculated the luminosity to 
0(v 4 ). Then extending Poisson's method, a more systematic method was developed by 
Sasaki []T2J (hereafter Paper I) and analytical expressions for the ingoing-wave Regge- 
Wheeler functions X 1 ^ were derived with the accuracy required to calculate gravitational 
wave forms and luminosity up to 0(v 8 ) beyond Newtonian. There it was also shown that 
as long as we are concerned with the radiation going out to infinity, the effect of the 
presence of the black hole horizon does not appear until we calculate to an extremely 
higher order, 0(v 18 ). Hence the expansion of X^ can be used in a situation in which 
there is a non-rotating compact star instead of a black hole. 

In this paper, using the result obtained in Paper I, we calculate analytical expressions 
for the gravitational wave forms and luminosity from a particle in circular orbit around a 
Schwarzschild black hole to 0(v 8 ) beyond Newtonian. This paper is organized as follows. 
In section 2, we show the general formulae and conventions used in this paper. In section 
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3, we first briefly review the result of Paper I and discuss the relation between the Regge- 
Wheeler and Teukolsky functions with emphasis on the consistency of the post-Newtonian 
orders in the conversion formula. In section 4, we present the analytical expressions for 
the gravitational wave forms and luminosity to 0(v s ) beyond Newtonian and compare 
the result with that of Tagoshi and Nakamura |n]] . We find an extremely good agreement 
between the two. In section 5, applying our result of the luminosity formula, we discuss the 
effect of the higher order post-Newtonian terms to the accumulated phase of gravitational 
waves from coalescing compact binaries. Section 6 is devoted to conclusion. Throughout 
the paper we use geometrized units, c—G—1. 



2 General formulation 

We consider the case when a test particle of mass /i travels a circular orbit around a 
Schwarzschild black hole of mass M >> [i. We mostly follow notation used by Poisson 
T3| , but for defmiteness, we recapitulate necessary formulae and definitions of symbols 



in this section. 

To calculate the gravitational wave forms and luminosity, we consider the inhomoge- 
neous Teukolsky equation [111, KB, 



dr \Adr v ; 



Rimu{r) = T em uj(r), (1) 



where 
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U(r) = ^[uV-4iw(r-3M)] - (£-l)(£ + 2), A — r(r — 2M), (2) 

and Ti muJ is the source term whose explicit form will be given in Eq.(^) below. 

We solve Eq. ([I]) by the Green function method. For this purpose, we need a homoge- 
neous solution Rfa(r) of Eq.(|l]) which satisfies the following boundary condition, 

oin (r) = f D £ul A 2 e~ iuir * for r* -> -oo , , , 

&A ) | r 3£out e iwr* + r -l£in e -iu,r* for r * _^ +QO ? KJ 

where r* = r + 2Mln(r/2M — 1). Then the outgoing-wave solution of Eq.(|I]) at infinity 
with the appropriate boundary condition at horizon is given by 



Rtmu(r -> oo) = . / drRZT lmu] {r)A- 2 

= r 3 e^ r *Z emuJ . (4) 

In the case of a circular orbit, the specific energy E and angular momentum L of the 
particle are given by 



E= (r -2M)/ v /r (r -3M) ) (5) 

and 



L = JMr /Jl-3M/r , (6) 
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where tq is the orbital radius. The angular frequency is given by Vt = (M/tq) 1//2 . Then 
Ttmu{r) is given by 



-2 b im {r - 2M) 2 5(r - r ) 

^h m 2ir Q \(r - 2M) 2 5\r - r ) - (r - 2M)(2 - zo;r )5(r - r ) 



+- 2 b, 



2^lm 



r J (r -2M)Y(r-r ) 



+{2icjr|j(r - 2M) - 2r (3r 2 - 8r M + 4M 2 )}5'(r - r ) 

+ [4r 2 - 8M 2 - u 2 r* - Qiu;r 2 (r - M)}d(r - r ) \d(u - mfi), (7) 

where <5(r) is the Dirac delta function and ' = d/dr. The coefficients s bi m are given by 



-ibem 

-2b im 



-[(£-!)£(£+ l)(£ + 2)] 1 / 2 Y em ( £,0) £r /(r - 2M) 



-f)(£ + 2)] 1 / 2 _ 1 F, m (^0)L/r , 



7T 



-l^m ( 2 j ) 



where s Yi m (6, (p) are the spin-weighted spherical harmonics. From Eqs.(f§) and (0), we see 
that Ze muJ takes the form, 

Zimuj = Z lm 5(u - mfi), (9) 



where 



7T 



2zo;r 2 J Bt 



- ^ m - 2i _!^ m ^1 + -wr*/ (r - 2M)^ 
+i _ 2 &£m^o(l - 2M/r )~ 2 f 1 - M/r + \iur 



+ 



i „ x bim - -2km ( 1 + -ujr /(r - 2M) 



+ -^2b em ro 2 RZ"(r )\. 



(10) 



We note the symmetry, Ze- m = (—i) Zg m , which may be seen from Eqs.(||), (|T0D and the 
symmetry of the spin weighted spherical harmonics, s Y^ m (7r/2, 0) = (— l)( s+ ^ s Y^ m (7r/2, 0). 
In terms of the amplitudes Ze m , the gravitational wave form at infinity is given by 



and the luminosity is given by 



dE 
~dt 



ira 



00 t 



;n) 



'12) 



1=2 m=l 
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where u = mVt. Thus the only remaining task is to calculate the ingoing-wave Teukolsky 



function R 1 ^. 



In stead of directly calculating R 1 ^ from the homogeneous Teukolsky equation, it is 
much easier to calculate the corresponidng Regge- Wheeler function Xg^ first and then 
transform it to Rf^. The homogeneous Regge- Wheeler equation takes the form [p|, 



dr 



- + <S-V{r) 



Xfjr) 



where 



V(r) 



r J 



0. 



6ikf 



It is known that this equation is obtained by transforming R^ as [|T7|, [18 



(13) 



(14) 



(15) 



Conversely, we can express X^ in terms of R^ by the inverse transformation formula as 



X, 



d 



r 2 ( d 



R 



c A \dr* y A \dr* 



luJ 



r 



2 ' 



(16) 



where Cq = (£ — 1)£(£ + + 2) — YIiMuj. Then we obtain the asymptotic forms of X t 
as 



m 



e ^ ^,.>C 



C^ 1 



-oo. 



(17) 



where A^,, and are respectively related to B^, B™' and defined in Eq.@ 

as 



75111 
75 OUt 



4a; 2 ^' 
-4u; 2 ^\ 



c 



16(l-2iMu;)(l-4iMw) 



3 Post-Newtonian expansion of the Teukolsky func- 
tion 



3.1 Method 

In Paper I, the post-Newtonian expansion of the ingoing-wave Regge- Wheeler functions 
X^ was formulated and they were calculated to 0(e 2 ) where e = 2Mu. In this subsection 
we briefly review the method. 

We rewrite the homogeneous Regge- Wheeler equation as 



1_ 

dz 



*2 



+ 1-1- 



+ 1] 



3e 



XI 



0. 



(19) 



5 



where z = wr, z* = z + eln(z — e) and we have suppressed the index uj since it is trivially 
absorbed in e and z. Note that for uo = mQ, e = 2mv 3 and z = mv at r = r . Hence the 
post-Newtonian expansion corresponds to expanding X™ with respect to e and evaluating 
X™ at z <C 1 as well as Af to required orders in e. Now setting 



X, 



in _ g — ieln(z — e) 



we find that Eq.(19) becomes 



d 2 _d_ / 

dz 2 z dz \ 



6 = ee lz — 
dz 



(20) 



(21) 



Thus expanding ^ with respect to e as 

oo 

&(*) = E^ B) ( 



(22) 



n=0 



we obtain the recursive equations, 



d 2d / 

dz 2 2; dz \ 



First for n = 0, we have 



+ 1)' 



d B) c*) 



ee 



dz 



1 d 

z 3 dz 



(23) 



(24) 



where j% and are the usual spherical Bessel functions. The boundary condition is that 



Q ' be regular at z = 
n > 1, we rewrite Eq.(B 



]. Hence /3 < -°- ) = and for convenience we set 
in the indefinite integral form, 

1 



ft 



(0) 



1. For 



rm I dzz 2 e tz ji 
—jt \ dzz 2 e~ tz ri(; 



(e lz z 2 ^ 



iz 2An-\), 



(25) 



If the above indefinite integrals can be explicitly performed, it is easy to obtain Q n ^ with 
a desired boundary condition. In Paper I, it is shown that this is indeed the case for n = 1 
and 2, and the boundary condition is that £^ be also regular at z = at least for n < 3. 
We then find for n = 1, 



(i) 



(£-!)(£ + 3) , / I 2 -4 2£-l 
2(£ + l)(2£ + l) Jm ~ ' :T777 ^ 



\2£(2£+l) £(£ 



Jt-i 



1-2 



+z 2 {n e j - j e n )j + ( 7 



fc=i 



fc fc + 1 



z 2 (n e j k - jen k )j k 



+ri£(Ci2z — 7 — In 2z) — Si2z + ijg lnz + oq ji 



(26) 



where Ci(x) = — dt cost/t and Si(x) = f£ dtsmt/t are the cosine and sine integral 



functions, respectively, and otp is an arbitrary integration constant which represents the 
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arbitrariness of the normalization of Xg n . We set a t — 0. As for n = 2, closed analytical 

(2) 

expressions of Q for £ = 2, 3 and 4 can be found in Paper I, §4.2. 

To obtain Xf 1 from ^ is straightforward. Decomposing the real and imaginary parts 
of £^ as 

d n) = /, (n) +^ n) , (27) 
we find the imaginary parts gf^ are given by 

y (0) = , gW= Je \nz, = + fP In z. (28) 

z 

These relations confirm that Xf 1 is real at least up to 0(e 2 ). In fact, inserting Eq.(|2TD 
into Eq.(^) and expanding the result with respect to e by assuming z 3> e, we find 

XT = e- i ^^h(j £ + e(fi 1) +igi 1) )+e 2 (fP+igf ) ) + ..) 
= z (j t + e/« + e 2 (/« + g?> In, - \j^zf) +... 

= z (j t + ejf + e 2 (/i 2 > + zf) + . . .) . (29) 

3.2 Calculation 

Once we have closed analytical expressions of <Q n , A™ can be readily obtained by gathering 
all the coefficients of leading terms proportional to e~ tz / z from the asymptotic forms of 
£e at z = oo. The explicit forms of A™ from £ = 2 to 4 which are correct up to e 2 order 
are given by 

2 2 P V 3 420 

vr 2 f 25 5 2 107. , n .1 \ , on , 

e- + eN 1 tt 2 H (7 + ln2)^ + ... , (30) 

2 1 18 24 210 W V / V ; 

{ . 13 . ,13 



3 2 P V 6 84 / 



x l-e- + e 



7T f 169 5 2 13 

<^ 1 IT 2 -I 



2 I 72 24 



7r 2 + -( T + m2)} +...), (31) 



Af = Ue-^ 2 ^exp(ze—-ze 2 ^-TT 
4 2 F V 60 13860 

vr 2 [22201 5 2 1571 



x 1 - e- + e 2 <^ + — 7T 2 + (7 + In 2U + . . . , (32) 

1 2 I 7200 24 6930 u ; r f 



For £ > 5, explicit forms of A^ 11 which are valid up to order O(e) are given by 



A in = i/+ 1 e - ie ( ln2e +7) 
1 2 



tt ief^l < 1 (*-!)(* + 3) 
. '2 2\ j feifc + j fcifc + *(*+!) 



(33) 
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As discussed in Paper I, and will be mentioned below, the above formulae are what we 
all need for Af in order to calculate the wave forms and luminosity to 0(v 8 ) beyond 
Newtonian. 

In addition to A 1 ™, we need the series expansion formulae of JQ n at z <C 1. To clarify 
their required orders of accuracy, we first consider the qualitative behavior of Xf 1 as 



0. We find 



x 



(in) 



(in) 



0(1) 
O(z) 



eO{z) 
eO(l) 



l {0(z) 



0(l)\nz} + e :i O(z- 1 ] 
0{z) hiz} + .. 



X 



(in) 



0{z 1 - 2 ) + eO(z £ - 3 ) + e 2 {0(/" 4 ) + 0(z 



lnz} + ...J (£>4). (34) 

(35) 



The above result tells us the accuracy of X™ needed to achieve 0(v 8 ) beyond Newtonian. 



For convenience, we set 



XI 



First, note that we must calculate xf ] to 0(z u ), Xf to 0(z 8 ), X\*> to 0(z 5 ) and X, 



oo 
n=0 

11\ 



(») 



(36) 



'(2) 



.(3) 

to 0(z 2 ). Then we see from Eq.(p6|) that we need the series expansions of Xf 1 at z = 
to the above required orders for £ = 2 to 0(e 3 ), for £ = 3 4, 5 and 6 to 0(e 2 ), for £ = 7 
and 8 to O^ 1 ), and for £ = 9 and 10 to 0(e°). 

Since we have closed analytical expressions of Xf ] and X £ (1) for all £, and xf ) for 
£ < 4, whose series expansion formulae can be trivially obtained, what we additionally 
need to do is to derive the series expansions of X^ 3 \ X5 and Xg 2 \ Taking account of 
the required accuracy mentioned above, we find that only the leading terms in their series 
expansions are necessary. First consider Xg and Xg . With the bondary condition that 



these be regular at z — 0, and using Eq.(p5|) with n = 2 and Eq. (p6|) , They can be easily 
calculated. We find 



X, 



(2) 



Z 



X, 



(2) 



7425 
4 



212355 



+ 0(z 6 



z 5 + 0(z 7 ), 



(37) 
(38) 



where possible ambiguities due to the choice of integration constants appear only at 
0(z 6 ) and 0(z 7 ), respectively, hence do not affect the coefficients of the leading terms, as 

(3) 

it should be so. In the same way, the leading behavior of at z — is found to be 



X, 



(3) 



319 
6300 



z 2 + 0(z 4 ), 



(39) 



where the aforementioned ambiguity appears at 0(z 4 ), hence again does not affect the 
result. In Appendix A, we show the series expansions of X^ up to the required orders. 

To calculate gravitational waves from the Teukolsky equation, we need to know Rf^iz) 
at z -C 1, which we calculate by using the conversion formula from Xf 1 to Rf, Eq.(|15|). 
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uRf = z{z 




In terms of xj, n \ the result is expressed as 

z 2 

z(z - e) 

= -z 2 {zXf y ) + z 2 (zXPY + 2zz 2 (zXP)' 
+ e (-2 z (zXP)" + z 2 (zX^)" - z 2 (zX^) 

-2 i z {zXf)' - i (zXf ] ) + 2iz 2 (zXf y ) 




+ e 2 (+(zXPY - 2 z (zXP)" + z 2 (zXP)" - z 2 (zXP) 

-i (zXf ] ) -2iz (zX^y + 2iz 2 {zXf ] ) 
+ e 3 (HzX?)" - 2 z (zXP)" + z 2 (zX®) 

-i (zXP )-2iz {zXf ] y + 2tz 2 {zXf ] y ) , (40) 



where ' = d/dz. 

Here a problem might arize, since what we have calculated are the series expansions 
of Xl n (z ) to 0(v 8 ) relative to the leading order, but what we actually need are those of 
R™. Hence we must examine if there is consistency of the post-Newtonian orders between 
the series expansions of X™ and R™. 

A straightforward way to check this consistency is to insert the truncated series ex- 
pressions of X^ accurate to 0(v 8 ) beyond Newtonian to Eq.(^), insert the resulting 
expression of R™ to the inverse transformation formula fll6|), and examine if the final 
result is the same as the starting expression of Xf 1 . We have verified this is indeed the 
case. 

As for B 1 / 1 , it should be now obvious from Eqs.(|TT|), (Q) an d ( |18D that the calculated 
Af have sufficient accuracy to obtain Bf 1 to the required orders. 

All these facts support our confidence that X} n and Rf 1 are equivalent also within the 
framework of our post-Newtonian expansion. 



4 Wave forms and luminosity to 0(v 8 ) 

Let us now calculate the gravitational wave forms and luminosity up to 0(v 8 ) beyond 
Newtonian. The task is straightforward but tedious. So we only show the key equations. 
The calculations have been performed with the help of the computer manipulation soft- 
ware Mathematica. We follow the notation of Poisson [|I^] to describe the post-Newtonian 
expansion of gravitational wave forms and luminosity. Namely, we express them as 

dE 32 /M 2 /M\ 5 » ~ 

if = tU) U (41) 



>.„> +h l~rn = -[^){^)Q m - (42) 



9 



4.1 Luminosity 

The procedure is straightforward. First we obtain Rf from Eq . fl40|) and Bf 1 from Eq . (|i"8|) . 
Then we insert them to Eq. ([To|) and set z = mv and e = 2mv 3 to obtain Z^ m . Finally, 
inserting them to Eq.([T2|) and expanding the results with respect to v, we obtain r\i m . 
The explicit forms of rj£ m are given by 

107 v 2 , 4784 v 4 428vrt> 5 19136 vrw 7 
772,2 = 1 ^ h4?iV + — — — h 



+v 6 



21 1323 21 1323 

/ 99210071 1712 7 16 tt 2 3424 In 2 1712 In?/ 

V 1091475 T05~ + ~~ 3 105 105 t 



27956920577 183184 7 1712vr 2 366368 In 2 183184 In v 
+v° I — — — — + — — '- — — + — — + 



81265275 2205 63 2205 2205 

v 2 17 v 4 ttv 5 2215 v 6 17 ttv 7 



772,1 36 504 + 18 254016 252 

8 / 15707221 1077 ^ 107 In 2 107 lnv\ 

+V \ 26 195400 945 + 27 945 945 J ' 

1215 v 2 1215 v 4 3645 tt 243729 v 6 3645 nv 7 

V3,3 = TTTT- H 777 h 



96 112 448 9856 56 

8 / 25037019729 _ 47385 7 3645 vr 2 _ 47385 In 2 

+V V 125565440 1568 + 224 1568 
47385 In 3 47385 In 



1568 1568 
5v 4 193 v 6 20nv 7 86111 



773,2 = — r"— + 



63 567 ' 63 ' 280665 ' 
v 2 v 4 ttv 5 437 v 6 nv 7 



^3,i = 7777 - 7777 + 77777: + 



8064 1512 4032 266112 756 



s , 1137077 137 tt 2 13 In 2 13 Inv 
~ ' ~ 50854003200 " 42336 + 6048 _ 42336 _ 42336 

10240 tt v 7 560069632 v 8 

//L1 " 567 6237 + 567 + 6243237 

729 v e 28431 v 8 

5?4,3 



1280 v 4 


151808v 6 


567 


6237 


729 v e 


28431 v 8 


4480 


24640 ' 


5v 4 


437 v 6 20 ttv 






3969 


43659 ' 3969 




101 v 8 


282240 


4656960' 



7199152 v* 

'1^ - 777^ _ 'io,:r„ 1 orw:n f 218513295' 



574,1 
575,5 
575,4 
575,3 



9765625 2568359375 w 8 



2433024 47443968 
4096 

13365 ' 

2187 v 6 150903v 8 



450560 2928640 
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V5,2 



4 v s 

40095' 

v 6 179 v s 

7/5,1 ~ 127733760 ~ 2490808320' 
26244 o 

^ = "3575 V ' 
131072 8 

V6A ~ 9555975* ' 

4 8 

' ' 5733585 



From the above results, we obtain dE/dt which is correct up to 0(v 8 ) as 

dE ( dE\ ( 1247 v 2 , 44711 v 4 8191 ttv 5 16285 nv 1 

+ 4 7T IT ~ 



dt \ dt J N \ 336 9072 672 504 

/6643739519 1712 7 16vr 2 3424 In 2 1712 In 



+v 6 



+v 



V 69854400 105 3 105 105 

323105549467 2325977 1369 tt 2 39931 In 2 



3178375200 4410 126 294 

47385 In 3 232597 In v\\ 



1568 4410 



(43) 



4.2 Wave forms 



Noting the form of the spherical harmonics, s Ye m (8,(p) = s Pe m (8)e imip , and that of A 
given by Eqs. (150) ~ fl3"3]), we find the wave form is in the form, 



htm-ihf m OC e «(la26+7) e -imfi(t-r-*) e ^ 



_ e im[2v i (-y+21n2+31nv)-n(t-r*)+<p] e 2imv i lnm 

This suggests that it is more convenient to introduce a new phase variable ip as 

V» = - r*) - y» - 2w 3 ( 7 + 2 In 2 + 3 In u). (45) 

An advantage of introducing the phase variable ip is that we directly see the post- 
Newtonian corrections at 0(v 3 ) to the phase of waves. These phase corrections at 0(v 3 ) 
delay the arrival time of the wave to the observer and express the tail effects. From 
Eqs. (|30|) ~ (|32|), we also see the tail effects at 0(t> 6 ) on the phase. 

In passing, we note that if one expanded A l f in terms of v, there would appear terms 
like (lnt>) 2 at 0(v 6 ). Then one would be surprized to see the cancellation of them in the 
calculation of 0(e 2 ). However in our approach, we see from the beginning that such terms 
do not exist. Together with the fact that there are no Inv terms at orders smaller than 
0(v 5 ) and no (lnv) 2 terms at all in Xf 1 , we readily foresee that there will be no (Inv) 2 
terms in dE/dt to 0(v s ), as we have explicitly seen in the previous subsection. 
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Using ip we obtain 

him + 
\i 1 



r uj 



( -2Pim(0) + (-1)' - 2 iVm(0)) [^ m e-^ + W^e' 



(46) 



u £m ^ '%-m ~ 

r uj v 

where W/n, is defined as 



-imif) 



W lm e 



imtp 



J im 



\xr Jmip 2u 3 (7+2 In 2+3 In 



(47) 



(48) 



Now from Eq. (|42|) and using the exploit forms of -2-Pfm(^) ; it is straightforward to obtain 
Ce~ m x - The explicit wave form for each £ and m are given in Appendix B. 

4.3 Comparison with the numerical result 

In this section we compare our analytical result for dE/dt with the numerical result 
calculated by Tagoshi and Nakamura |TT] . 



The numerical values of the coefficents in the luminosity formula calculated in this 
paper are given in Table 1(a). Comparing them with the coefficients given in their paper, 







v 4 


-4.928461199295258 


v 4 


-4.928461199294533 




V 5 


-38.29283545329089 


V 5 


+38.29283545469344 




V 6 


+ 115.73172132 


V 6 


+115.7317166756113 




v 6 Inv 


-16.304761151 


v 6 lnt> 


-16.3047619047619 




V 7 


-101.509987 


V 7 


-101.5095959597416 




V 8 


-117.787 


V 8 


-117.5043907226773 




v 8 \nv 


+52.6901 


v 8 Inv 


+52.74308390022676 




Q 

V 


+700.45 






v 9 \nv 


-209.78 




(a) 




(b) 



Table 1: (a) The numerical values of the analytically calculated coefficients of the post- 
Newtonian expansion of dE/dt given in Eq. (|43|) , and (b) those numerically calculated by 
means of least square fitting, where the same numerical data as in Tagoshi and Nakamura 
llj are used, but additionally including \nv terms at 0(v 9 ) and 0(v w ). 



we see that their coefficient of v 8 differs from the true value about ~ 20%, while all 
the other coefficients are in very good agreement with each other. This difference is 
within the error estimated in their paper. However, we find that in fact their numerical 
data turn out to have much better accuracy. The reason is as follows. In Ref.|ll|, the 
coefficients are calculated by least square fitting, but taking account of In v terms only at 
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0(v 6 ) and 0(v 8 ). However, since they calculated the gravitational waves to £ = 6, their 
data correctly contain contributions from the orders up to 0{v 9 ) of the post-Newtonian 
expansion and some constributions from yet higher orders, provided the data have enough 
accuracy. Hence we have recalculated the coefficients using their numerical data by least 
square fitting, including v 9 lnv and t> 10 lnt> terms. The result are given in Table 1(b). 
We find that even the coefficient of v 8 agrees with the analytical value within 1%. This 
suggests that their data do actually give the coefficients to 0(v 9 ) of the post-Newtonian 
expansion with high accuracy, hence the values of the coefficients of v 9 and v 9 In v given 
in Table 1(b) are expected to be good approximations to the true values. 



5 An estimate of orbital phase of coalescing binaries 

In this section, based on our results, we discuss the accuracy of the post-Newtonian expan- 
sion to construct the template wave forms from inspiraling compact binaries. Although 
our results are valid only for the test particle limit, ji/M <C 1, we ignore this fact in 
the following. Since the effect of no n- vanishing fi/M would only increase errors in the 
estimate, our estimate below may be regarded as an optimistic one. 

The total cycle N of gravitational waves from an inspiraling compact binary during 
sweep through, say, the LIGO band is 

where / is the orbital frequency, £«(£/) and rj(rj) are the initial (final) time and the 
orbital separation of the binary, respectively, and we have assumed quasi-periodicity of 
the inspiral orbit. Then expanding both the numerator and denominator with respect to 
v, N is expressed as 

5 M r*s . ET=ohx k 



N = 



327T /i Jxi x 6 EfcLo a fc xfc ' 

where x = (M/r) 1 ^ 2 and the series forms in the denominator and numerator sympoli- 
cally represent the post-Newtonian corrections to the dE/dr and dE/dt, respectively, i.e, 
Y,b k x k = (dE/dr)/(dE/dr) N and £ a k x k = (dE/dt)/{dE/dt) N . 

To examine the accuracy of the post-Newtonian expansion, we introduce which 
is defined by 

Nin) = J_M. [ Xf dx (51) 
32vr fj, J Xl x 6 ELo a kX k ' 

To be specific, we assume the detectable frequency band to be from 10 Hz to 1000 Hz 
and the quasi-periodic inspiral stage ends at r = 6M. So, rj is the radius at which 
fin) = 10 Hz, and r f is the one at which f(r f ) = 1000 Hz if r f > 6M and r = 6M 
otherwise. 

For (1.4M , 1.4M ) neutron star binary, we have = 175M and rj = 8M. Then 
we obtain = 16220.13, = 16211.03, = 16211.93 and = 16211.91. For 
(1.4M , 1OM ) binary, n = 68M and r f = 6M, and we obtain iV^ = 3484.83, iV"( 6 ) = 
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3466.56, = 3468.60 and = 3468.27. For (1OM , 10M o ) binary, n = 47M and 
T f = 6M, and we obtain = 581.39, = 574.49, = 575.29 and = 575.15. 

From these results we see that for (1.4M , 1.4M ) binaries, the post-Newtonian ex- 
pansion to 0(v 7 ) will be necessary to accurately predict the total cycle. However for 
(1.4M , 1OM ) and (1OM , 1OM ) binaries, the convergence seems to be slow. The error 
in the total cycle seems marginal, AiV < 0.5, even including up to 0(v 8 ). Further, if we 
regard the value of the coefficient of v 9 in Table 1(b) to be approximately correct, about 
700, the contribution from 0(v 9 ) may be comparable to 0(v 8 ). For these binaries, more 
detailed analyses are needed. 



6 Conclusion 

Based upon formulae developed in Paper I, we have analytically perfomed the post- 
Newtonian expansion of the gravitational waves from a test particle in circular orbit 
around a Schwarzschild black hole. 

We have calculated both the gravitational wave forms and luminosity to order v 8 
beyond Newtonian. In particular, we have given the exact analytical values of the coef- 
ficients of lnt> terms at 0(v 6 ) and 0(v 8 ) in the luminosity. The existence of such terms 
was found previously in a numerical analysis by Tagoshi and Nakamura [pT[| . We have 
found that their numerical results are in very good agreement with our analytical ones. 

We have also estimated the accuracy of the post-Newtonian expansion to predict the 
total cycle of coalescing binaries. For (1.4M , 1.4M ) binaries, we have found that the 
post-Newtonian expansion up to 0(v 7 ) will be sufficient to construct accurate theoretical 
templates. However for (1.4M ,1OM ) and (1OM , 1OM ) binaries, accuracy to order 
higher than v 8 seems to be necessary. 

Although valid only in the test particle limit, our analytical results should be repro- 
duced in the conventional post-Newtonian calculations which are not restricted to the 
test particle case. Hence our results give a useful guideline for the future researches in 
the gravitational wave physics of coalescing compact binaries. 
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A The expansion of X 



in 



In this Appendix, we give the expansion forms of which are required to calculate 
Xf 1 up to order v 8 . 



x (o) = ^__J^ + _ z 9 + : 1 



15 210 7560 498960 51891840 



14 



-13 z A 



+ 



630 810 
3 / 26743 

Z \ 110250 ~ 

319 z 2 



53 z 8 

~ 1782000 
107 \nz 

3150 



+ 



227 z 



10 



6300 

-v4 



+ z a 



AO 



567567000 
140953 



+ 



107 Inz^ 



9261000 44100 



105 



390 

v 5 



+ 



+ 



83160 
221 z 1 



6486480 



126 630 2494800 
76369 2 s 2327663 z 5 



1852200 1100206800 

^5 SI <A 



13z 3 \nz 13z 5 lnz 
+ 



4410 



79380 



945 20790 



+ 



— z 

630 

z 3 

1764 

z 6 

10395 
-z b 

4950 

-v4 



Z 



+ 



1081080 

79 z 8 



97297200 



9900 13899600 

(942578261 - 43548120 In z) 

181534122000 

z 8 z 10 



270270 

41 z 7 

8108100 



16216200 



7425 

z 7 

135135 
-8z 6 

405405 

4z 5 

212355 



2027025 

-z 7 

630630 

9 

Z 

34459425 " 

-z 8 

9189180 

z 10 

654729075 



z 



4054050 

z 8 



+ 



z 



11 



275675400 



5675670 



,10 



68918850 



,ii 



1309458150 



13749310575 



15 



B Wave forms to 0(v 8 ) 

In this Appendix, we give the gravitational wave forms for all the relevant £ and m that 
contribute up to 0(v 8 ). 
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_ (cos(fl) + 7cos(3fl)) / , 437 v 4 sin(2^) 

C 4 , 2 - ^ \v sm(2^ — 

148 \ . \ 1038039v 6 sin(2^) 



+v 5 (cos(2 ^) (— - 4 ln 2 j + 2 tt sin(2 + 

7 / /rt / 32338 874 1n2\ 437tt sin(2^)\ 

+ .„ 8 ( cos(2 . W (^£_ 8lln ^ ' 54548715967 



200200 



2497294800 



125687 2tt 2 111616 ln2 x2 12568 \nv\ . " 

+ + 8 In 2 2 s in 2 ib) 

3465 3 3465 V ; 3465 J V ^, 

_ (-2 sin(2g) + 7sin(4fl)) / 101 v 5 cosgj 

C4 ' 1_ 560 r C ° SW 33 



23 



42982v 7 cos(^) «( 149sin(^)" 

+ 15015 + " 30~^ 

8 / -101 7T cos(^) 15049 sin(V>) \ \ 
+V [ 33 + 990 J J 



_ 625 (-2 sin(2fl) + sin(4fl)) / 263 v 5 cos(5^) 

C 5 ,5 - 7^ \v cos(5^j - 

569 m i A • /c 9185 cos(5^) 



819 



/ / 569 \ \ 

+v 6 f 5 7T cos(5V) + ( — — + 10 In 5 J sin(5^)J + 

w j-1315^) + (mm 2^5 j 

-2 (14 cos(0) -9 cos(3 0) -5 cos(5 0)) / 4 , . 4451 v 6 sin (4 tb) 

Cm = — 4^ — {«' M**) 

W (cos(4^) (f - 8 .„4) +4, Si „(4 ,)) + 10715 2 "; 8 7' 4 ^ ) 



& = -27(2^(2^-9^(4 0)) / 3 cos(3v , ( _ 69 v" cos(3*) 



1280 V 13 

569 ci o\ ■ (o ,\\ 12463 w 7 cos(3^) 



1365 



/ / o69 \ 

+v 6 i 3 7T cos(3 ^) + ( -— + 6 In 3 J sin(3 ^) 

o /-207tt cos(3^) /39261 414 ln3\ , ,.\\ 

+v { vT^ + ( -455- is-j Sln(3 ^jj 

(14 cos(fl) + 3 cos(3#) + 15 cos(5#)) / 4 . , rt , N 3911 v 6 sin(2 ^) 

(, x , = — — — — v sm(2 ib) —- - 

S5 ' 2 540 V 910 

7 / /« ,x /163 , , rt \ rt . , rt , x \ 63439 v 8 sin(2^)\ 
+v 7 (cos(2 ^) (— - 4 In 2 J + 2 tt sm(2 V) J + — 1 

(2 sin(2#) + 3 sin(4 0)) / „ , , N 179 v 5 cos(^) 

C= , = 1 — 1 — — i» cos(^) — 

S5>1 5760 V 39 

fi / , 569 sin(^)\ 5023 v 7 cos(^) 

+v 6 ir cos(ib) — H — 

V K J 105 J 585 

8 /-179vr cos(^) 101851 sin(^)\\ 
+V { 39 + 4095 J J 

81 cos(0) sm 4 (9) / , . 113 sin(6^) 

Ce ' 6 = 10 K Sm(6 ^ 14^ 

487 \ . \ 1372317 v 8 sin(6V)\ 



+v 7 (^cos(6 ^) (— - 12 In 6^ + 6 tt sin(6 + 



73304 



-3125 cos(#) (1 + 3cos 2 (#)) sin 3 (#) / - . 149v 7 cos(5V) 

C 6 *5 = \mr — — V COS ( 5 ^) ~ 



2016 V v ' 24 
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/ / 1219 \ 

+v 8 i 5 n cos(5^) + f — + 10 ln5J sin(5^) 



-64 cos (9 1 + 11 cos (29 sin 2 (9 ) 4 93 v 6 sin 4 ^ 
, fi 4 = 1 — — v smUib) y -— L 

7 / n /487 \ A . IA ,,\ 3261767v 8 sin(4^)~ 
+v 7 cos(4 V) — - 8 In 4 + 4 7r sm(4 if>) + 



,3 



21 ) v r V 219912 

243 cos(#) sin(0) (53 + 20 cos(2 0) + 55 cos(4 9)) / 5 

gg560 (u cos( 

133 cos(3^) 8 / / 1219 



24 (3ttcos(3^) + (-— + 6 1n3) sin(3^)) 

_ cos(0) (47- 84 cos(2 0) + 165 cos(4 0)) , 4 

oo 7fin r S1I H 2 ^ 



23760 

81 v 6 sin(2^) 7 / , rt /487 



+ v 7 ^cos(2^) 
+2tt sin(2^)) + 

^6,1 



14 V V 42 

14482483 v 8 sin(2^)' 



4 ln2 



1099560 J 
cos(0) sin(0) (41 - 12 cos(2 0) + 99 cos(4 0)) 



266112 

' , 125v 7 cos(vV) 8 / , 1219 sin(^)' 

V 5 COS(^) + V S 7T COS(^) 



117649 cos (0 sin 5 (9 ) = ,„ 319 w 7 cos 7^ 

C 7 7 = — — cos(7^) ^— ^ 

S7 ' 7 11520 V 34 

+v 8 (Vtt cos(7^) + + 14 ln7 ) sin(7^))) 

_ 243 cos(fl)sin 4 (fl) (13 + 7 cos(2fl)) 
Ct ' 6 ~ 560 

/ fi , 1787 v 8 sin(6^)\ 

-78125 cos(fl)sin 3 (fl) (3 + 13 cos(2fl)) , 
C 7l5 - [y cos(5^ 

271 v 7 cos(5V) , .. 8 f c _ _, C .,. N , f 7699 
34 

-4cos(#) (113 + 116 cos(2#) + 91 cos(4 #)) sin 2 (fl) 
C? ' 4 ~ 1365 
f e . 14543v 8 sin(4^)\ 

V v 2142 ; 

-729 cos(#) (167 + 44 cos(2#) +429 cos(4#)) sin(0) 

10250240 



/ / 7699 \ 

+ v 8 [5 7v cos(5 tp) + f — 2^2" + 10 ln 5 ) sin ( 5 ^) 
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239 v 7 cos(3V) 



v cos(3-?/>) — 



34 

7699 



+V 8 (3 7T C0S(3 V) + (—420" + 6 111 3 ) Sill ( 3 ^)) 

cos(#) (-338 + 1351 cos(2 9) - 990 cos(4 6) + 1001 cos(6 6)) 



7,2 768768 
v 6 sin(2^) 



13619 v 8 sin(2^)\ 
2142 J 
cos(6) sin(0) (109 + 132 cos(2#) + 143 cos(4fl)) 
18450432 

5 223v 7 cos(^) 8 / 7699 sin(^)' 



IT COS(V') h l> { 7T COS(V') 



34 V 1260 

ff 8 = -««Wco.ffl.l.«ffl / a sin(g ^ 3653 S in(8 > 



315 V 342 

!? 

207360" 



= fn^ 7 COS ( 7 ^ (2 + cos(20)) sin 5 (£) 



_ 729 cos(fl) (1 + 3 cos(2fl)) sin 4 (fl) 

^ 8 ' 6 ~~ 560 

/ 6 353 1; 8 sin(6^ 
I u° sm(6^) — I 

— 781 2 s ) v 7 

Q 5 = n61216 cos(5^) cos(#) (11 + 14 cos(2#) + 7 cos(4#)) sin 3 (#) 

-4 cos(#) (49 + 52 cos(2 0) + 91 cos(4 6)) sin 2 (0) 
^ 8 ' 4 ~ 4095 
f 6 . n 2837 1; 8 sin(4^)\ 

^ Sm ( 4 ^ 342 J 

243 ?; 7 

Cs ' 3 = 3727360 008(3 ^ COS( ^ } (24 + 141 COS(2 ^ } + 91 C ° S(6 ^ Sin( ^ 
_ cos(0) (-274 + 583 cos(2 6) - 286 cos(4 6) + 1001 cos(6 9)) 

Q,2 — ~~ 



1 1531520 
2633 v 8 sin(2^V 

v & sin(2^) - 



342 



t; 7 



C & = 94887936 °° S( ^ ( " 8 + 121 C ° s(2 ^ + 143 C °^ 9)) ^ 

= -4782969 v 7 ^ ^ 7 ,^ 

9,9 286720 
— 8192 ?; 8 

C^s = 14175 ^(8^) cos(0) (19 + 9 cos(2 0)) sin 6 (0) 
403^3607 v; 7 

Cs ' 7 = 112803840 C ° S(7 ^ C ° S( ^ (7+ 17 C0S(2 ^ )} Sin5( ^ 

C 9 * 6 = 729 4 8 7 g n (6 ^ ) «»(*) (97+ 132 cos(2 0) + 51 cos(4 0)) sin 4 (0) 
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— 39062 c >7; 7 

Q 5 = 157Q25376 cos(5^) cos(tf) (59 + 80 cos(2£) + 85 cos(4£)) sin 3 (£) 

-4v* ° 
34425 



^4 = TT^sin(4^) cos(0) (142 + 343 cos(2£) 



+ 130 cos(4#) + 153 cos(6#)) sm 2 (9) 
243 v 7 

C^s = 7241728Q cos(3^) cos(^) (18 + 195 cos(2£) 
+78 cos(4#) + 221 cos(6#)) sin(0) 

Cs * 2 = H4566400 Sin(2 ^ C ° S{6) (1279_ 1480 cos ( 2 ^ + 
2236 cos(4#) - 1976 cos(6#) + 1989 cos(8#)) 

C 9 *i = n 73 ^QQaeo cos ^ cos ^ ( 166 + 403 cos ( 2 B ) + 234 cos ( 4 *) 

+221 cos(6#)) sin(0) 

, x 390625v 8 sin(lO^) . 8 

Cio.io = ' cos W sm W 

Cfo,s = ""^" (8r) «»(*) (9 + 19 cos(2 0)) sin^) 

Cio,6 = 218 1 Xo4o ( o V;) cos ^ ( 2449 + 3604 cos ^ + 2907 cos ( 4 ^)) sin ^) 

CiV = " 1 45^ 5 n 2 ( 5^ ) cos W ( 446 + 1319 cos ( 2 ^) 
+850 cos(4fl) + 969 cos(6#)) sin 2 (#) 

A* = ^ sm ( 2 ^) cos(^) (3007 - 5720 cos(2 0) 
S10 ' 2 17414092800 v ; v v ; 

+8268 cos(4#) - 1768 cos(6#) + 12597 cos(8#)) 
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